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ABSTRACT 

Molecular clouds often exhibit filamentary or sheet-like shapes. We compute 
the free-fall time {ts) for finite, uniform, self-gravitating circular sheets and fila- 
mentary clouds of small but finite thickness, so that their volume density p can 
still be defined. We find that, for thin sheets, the free-fall time is larger than 
that of a uniform sphere with the same volume density by a factor proportional 
to VA, where the aspect ratio A is given hy A — R/h, R being the sheet's ra- 
dius and h is its thickness. For filamentary clouds, the aspect ratio is defined as 
A = L/TZ, where L is the filament's half length and TZ is its (small) radius, and 
the modification factor is a more complicated, although in the limit of large A it 
again reduces to nearly y/A. We propose that our result for filamentary shapes 
naturally explains the ubiquitous configuration of clumps fed by filaments ob- 
served in the densest structures of molecular clouds. Also, the longer free-fall 
times for non-spherical geometries in general may contribute towards partially 
alleviating the "star-formation conundrum" , namely, that the star formation rate 
in the Galaxy appears to be proceeding in a timescale much larger than the total 
molecular mass in the Galaxy divided by its typical free-fall time. If molecular 
clouds are in general formed by thin sheets and long filaments, then their relevant 
free-fall time may have been systematically underestimated, possibly by factors 
of up to one order of magnitude. 

Subject headings: ISM: clouds — ISM: structure 
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Introduction 



The so-called free-fall time is one of the most important quantities in astrophysics. 
For a spherical object of mass M and radius i?, this timescale is given by (see, e.g., 



Binney fc Tremaine 



19871 ) 



37r 



SGM V 32Gp' 

where in the second equality we have introduced the volume density defined by 

3M 



p{M,R) 



(2) 



The timescale rg has the interesting property that it depends on the object's size and 
mass only through a combination that is proportional to its volume density, p. That is, 
once p is specified, Tg is independent of the object's mass (or size), implying that, in a 
collapsing uniform-density sphere, all spherical shells reach the center at the same time. 
This is equivalent to the well-known property that, for spherically-symmetric perturbations 
of a uniform medium, the growth rate increases with increasing w avelength, and thus th e 



fastest mode of collapse is an overall contraction of the medium (ITohline 



1980 



Larson 



19851). 



However, this independence of rg from the actual physical dimensions of an object of 
fixed volume density is only valid when the object's extension R is comparable in all three 
spatial dimensions (a "3D object"), because only in this case is the volume density of the 
object given by eq. Instead, for nearly sheet-like ("2D") or filamentary ("ID") shapes, 
the volume density p is not proportional M/L^, where L here generically denotes the 
object's largest dimension. For these morphologies, p oc M/{£L^) or pM/^M"^), respectively, 
where i denotes the fixed, small dimension(s) of the object. This is relevant for interstellar 
structures, since they are of ten observed to have sheet-like or filanientary, rather than 



spherical, morphologies (e.g. 



Bally et al 



1989 



de Geus et al. 



1990 



Heiles fc Troland 



2003 
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Mversl l2009l : iMolinari et al.ll2010l : lAndre et al.ll2010[ ). This suggests that the free-fall 
timescale may actually depend on the size of a non-spherical object in addition to depending 
its volume density. 



The gr avitational stability o: 



works (e.g., 



Ledoux 



1951 



Larson 



non-spherical str uctures has been considered in earlier 



1985 



Curry 



20001 ). but mostly considering infinite media 



and without discussing collapse timescales. Finite-size non-spherical structures, and their 
corresponding collaps e time s have only recently begun to be considered. In particular, 

I I I 

Burkert fc Hartmannl fl2004i) computed an approximation to the free-fall time rg^2D for 
finite-sized, infinitely thin circular sheets of radius R, given by 



'TS,2D ~ TfT,BH 



R 



(3) 



where S is the surface density (with units of mass per unit area) of the sheet. It is 
noteworthy that, in this case, the free-fall time exhibits a dependence on its size in addition 
to depending on the (column) density. Indeed, numeric al simulations of the collapse of 



Vazquez-Semadeni et al. 



(120071 ) 



large sheet-like clouds containing many Jeans masses by 
exhibited collapse timescales significantly larger than their corresponding three-dimensional 
free-fall time, as given by eq. ([I]). Motivated by these realizations, in this paper we comput e 



i n det ail the free-fall time for sheet-like and filamentary structures. Note that 



Pon et al. 



( I2OIII ) have recently investigated the free-fall timescales for sheet-like and filamentary 



geometries, although they have focused on whether small-scale perturbations within such 
structures have sufficient time to collapse before the whole structure does so. Here we 
concentrate on a different question: whether the collapse timescale for these geometries is 



longer, and 
and that by 



jy what arnount , compared to their spherical counterparts. Thus, our study 



Pon et al. 



(120111 ) can be considered as being complementary to each other. 
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2. The spherical case 

The sta ndard calculation of the fr ee-fall time for a uniform-density sphere proceeds as 



follows (e.g. 



Binney &: Tremaine 



19871 ). Consider a uniform sphere of radius R, mass M, 
and volume density p = P3d{M,R) = 3M/AitR^, which at time t = starts to contract 
under the action of its self-gravity exclusively. The subindex "3D" denotes the function 
that is used to calculate a quantity for a spherical (or "3D") geometry, which we distinguish 
from the physical value of the quantity itself, denoted without subscript. At a certain later 
time t > 0, when the sphere has radius r < R, the velocity of a point at its periphery is 



* -J2Gm("4). (4) 



dt \j \r R 

Introducing the non-dimensional variable x = r /R, this is equivalent to 



dx 
'dt 



2GM 




i?3 





(5) 



It is then clear that, because for a spherically symmetric uniform object the mass M 
increases as R^, expression depends on the mass and size of the sphere only through the 
volume density they imply, p = p-^D. Equation can then be integrated to yield 

™ - r.MM, R) ^ (^) /' (^) ^ (6) 

evidencing the independence of rg on the cloud's size or mass for a given volume density. 



3. Circular sheet-like cloud 



Let us now consider the case of an infinitely thin circular sheet of mass M, initial 
ra dius R, and uniform surface density S = M/tiR^. Here we follow closely the analysis 



by 



Burkert &: HartmannI (120041 . hereafter BH04). They calculated the radial acceleration 



experienced by the periphery of the sheet when it has shrunk to radius r < R, under 
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the assumption that the surface density S remains constant. This assumption was made 
because they found that the maximum acceleration occurs at the periphery of the sheet, 
and thus a dense contracting ring forms in the periphery, while the inner surface density 
remains essentially unchang ed. This behavior has been obse rved in numerical simulations 



by BH04 themselves and by 
found 



Vazquez-Semadeni et al. 



4GS 



R 



K 



E 



mm . With this assumption, BH04 

(7) 



.RJ \R 

where K and E are respectively the first and second elliptic integrals. Upon a series 
expansion, eq. ([7]) becomes 



vrGS 



r 3/r\3 45 /r\5 
^ ^ 8 \r) ^ 192 \r) ^ 



Note that the acceleration diverges at r = i?, but BH04 noted that this problem is 
eliminated when a finite sheet thickness is considered. 

BH04 retained only the linear term in eq. ([H]) and computed the free-fall time by noting 
that aj. = dv/dt = 1/2 dv^ jv dt = 1/2 dv^ jdr^ where v is the instantaneous radial velocity. 
They thus found the instantaneous velocity as a function of the instantaneous radius of the 
sheet, given by 



27rGS 



r'dr' 



ttGE 



R Jr R 
implying that the radial velocity at the end of the collapse is 



- r^) 



v{r = 0) = VnGJ^R. 



(9) 



(10) 



Taking this as a representative velocity for the entire collapse, BH04 then found a lower 
limit to the time for the sheet to shrink from r = R to r = 0, which we label rg^BH, given by 

R 



■^ff,BH 



v(r = 0) 



R 



as anticipated in eq. ([3]). 
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At this point we can extend the calculations by BH04 in two ways. First, we consider 
clouds with small but finite thickness, so that it is possible to define a volume density 
within them. We refer to this as a "quasi-2D" geometry. Specifically, we rewrite the surface 
density assuming a small but finite thickness h such that S = ph, where p is the volume 
density. However, in this case the volume density is not given by eq. ([2]), but rather by 
P = P2d{M,R) = M/irR'^h. Here, the subindex "2D" now denotes the relevant function to 
compute the quantity for a quasi-2D structure. 

Second, instead of taking the final velocity as representative of the entire collapse, 
we can obtain a more accurate expression by integrating eq. ([H]) from to r to write an 
expresion for v{r), obtaining 



Vi{r) = ^^iR'-r^), (12) 
where A = R/h is the aspect ratio. After a second integration we obtain 



n{t)=R sinl^-^^t), (13) 

where the subindex '1' denotes the assumption that the density of the sheet internal to its 
periphery remains constant. Setting ri = 0, we obtain the corresponding free fall time as 



rfr,i = W^ = VT"^'3D, (14) 

where the second equality compares with the free-fall time that would be obtained for a 
spherical structure with the same volume density, explicitly exhibiting the extra factor 
oc y/A. 

Let us now assume that, instead of remaining constant, the sheet's density 
increases during the contraction so as to maintain the sheet's mass constant; that is, 
T,{r{t)) = M/TTr{ty = p{r{t))h, with M = est. Note that, in reality, the sheet's mass 
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does remain constant during the collapse, but it is not distributed uniformly on the sheet. 
Instead, the mass external to r is piled up at the periphery. So, a point at the periphery 
sees this mass at the farthest possible distance within the sheet, rather than seeing it 
uniformily distributed over the sheet. Thus, the constant-mass assumption overestimates 
the gravitational pull of the sheet on this point. On the other hand, the const ant- density 
assumption neglects the mass at the periphery altogether. Thus, the two assumptions 
should bracket the real situation (still within the linear approximation). 

Under this assuption we can integrate eq. ([HD from i? to r to obtain a velocity function 

as 

.M^^2.Gp.«.n(f)^^M^>„(f), (15) 

where the subindex '2' denotes the case of a constant-mass assumption. Upon a second 
integration, we obtain 

r2{t) = ^ (16) 



exp 



where erf ^ is the inverse error function. This gives a free-fall time of 



rff,2=W^ = V^^ff^3D, (17) 

which is ~ 25% shorter than Tg^i because the gravitational acceleration is always larger, as 
the surface density increases monotonically as the sheet contracts. The estimates for the 
free-fall time in the two cases differ by factors of order unity at most. More importantly, 
both of them have the same dependence on the aspect ratio as ~ \/A. This dependence 
is fundamentally different from that of the spherical ("3D") case, given by eq. ([6]), as 
discussed in Sec.[Tl This is illustrated in Fig^{Top panels), which shows the infall velocities 
as a function of the instantaneous radius, corresponding to eqs. (fT^ and (IT3]) . taking 
R = p = G = 1, for a range of values of A. Also shown is the infall velocity for the spherical 
case. We see that the latter is always larger than either A values in both cases. This result 
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is clearly a consequence of the much smaller mass contained in a fiat, essentially "2D", 
configuration than in a spherical (3D) one of the same volume density. 

The bottom panels of Fig. [T]show the sheet's radius as a function of time for the two 
density cases. As expected, case 1, in which the cloud's mass is not conserved, but rather 
decreases in time, undergoes a slightly slower collapse than case 2. Nevertheless, we see 
that the difference in the collapse time between the two cases is around ~ 25%, and so both 
are qualitatively similar. 



4. Filamentary cloud 

We now consider the case of a uniform, cylindrical cloud of total length 2L, radius 
TZ L, and volume density p = pid{M,L) = M/tcTZ'^L, where now the subindex 'ID' 
denotes the function appropriate for calculating a physical quantity in the quasi- ID case. 
Note that here the filament's radius TZ is the small, fixed dimension, and the size variable is 
L. 



Again, we start with the case in which the filament's density remains constant. In this 
case, the acceleration towards the filament's center at a distance / from the center is given 
by (BH04) 



1 dv"^ ^ 



2/ + 7^ - ^fW^^ 



Integrating this equation from L to I, we obtain the radial velocity after the filament has 
contracted from L to I: 



v{l) = ^/AnGp 



I 



7^2 



{L-l){L + l + n)- -V^2:f4Z2 + ,/R2 + AP - — In 

^ ^ 4r 



(19) 



1/2 



If we now define the nondimensional parameters A = L/TZ and x = l/L, we can rewrite eq. 



1.0 0.8 0.6 0.4 0.2 0.0 

r/R 



1.0 0.8 0.6 0.4 0.2 0.0 

r/R 




00.00 



ff,3D 




00.00 



ff,3D 



Fig. 1. — Top: Evolution of the radial velocity of a sheet's periphery (shown as a function of 
its instantaneous radius, normalized to its initial value R) as it collapses, compared to the 
corresponding velocity for a spherical cloud (thick hne). Note that the radius decreases to 
the right. Bottom: Normalized position of the sheet's radius as a function of time, in units 
of the spherical free-fall time. Left: Case 1. Right: Case 2. The different lines show different 
values of A. 
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(IIHD as 



dx 



dt 2ts, 



3D 



+4 + 2 



^ + 4x2 + 2x 
(20) 



This equation again shows that the velocity depends on the "aspect ratio" A = L/TZ, 
besides the standard dependence on density given by rg 30- 

Let us now drop the constant-density assumption and assume instead that the volume 
density increases as the filament contracts as p{l) = Pid(M, /) = M/ttTZ'^1. In this case, we 
can perform the same calculation for the velocity evolution, to obtain 

1/2 

2(1 -x) - a/^ + 4 + a/^ +4x2 - ^In ' 'l = . (21) 



dx n 3 
'di ~ 2rff,3D V 2 



+ 4 + A/ + 4a;2 - T In 
\ A^ A 



4 + i 



1/2 



In Fig. [2](top panels) we show the contraction velocity as a function of / for a range 
of values of A, for both the const ant- density and constant-mass cases. In both, the final 
velocity reaches smaller values as larger values of A are considered. The lower panels in Fig. 
[2] show the position of the filament's edge as a function of time, obtained from numerical 
integration of eqs. fl20|) and fl2T|) . The results are normalized to the spherical free-fall time. 

Expressions fl2U]) and fl^ are quite complicated, and thus it is worthwhile to consider 
the limiting case of large A, for which an analytical solution can be found. We obtain 



2 8A 

TT V O 



and 



^5,2 




(22) 



(23) 



from which we see that, at large A, the free-fall time for filaments also exhibits an additional 
factor oc ^/A (Fig. [3]), similarly to the case for sheets. 



1.0 0.8 0.6 0.4 0.2 0.0 

X, r/R 



1.0 0.8 0.6 0.4 0.2 0.0 

X, r/R 




Fig. 2. — Top: Evolution of the edge velocity for a filamentary cloud, compared to the 
corresponding velocity for a spherical cloud (thick line). Bottom: Normalized position of 
the filament's edge as a function of time, normalized to the spherical free-fall time estimate. 
Left: Case 1, with p — est. Right: Case 2, taking p{l) — M/nTZ'^l. The different lines show- 
different values of A. 
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5. Discussion and applications 

5.1. Implications for filamentary structure in molecular clouds 

Our results have a number of implications for various aspects of interstellar structure 
and star formation. First, the fact that both the sheet-like and filamentary geometries have 
free-fall timescales larger than those of their three-dimensional (spherical) counterpart with 
the same volume density implies that any sheet or filament containing many Jeans masses 
will collapse later than any approximately-3D clump within it containing one Jeans mass. 
That is, the 3D clump will collapse on roughly one spherical free-fall time, while the sheet 
or filament around it will terminate its own collapse after a time that is longer by a factor 
This naturally explains the common 



within, and accreting 



2010 



Pillai et al 



rom, filaments (e.g. 



y observed morphology o 


' clumps immersed 


Mvers 


2009; 


Purcell et al. 


2009; 


Schneider et al. 



201ll ). because the clumps evolve on a shorter timescale than their parent 



filaments. 



Specifically, our results can be applied to t he filaments observed in the Aquila rift 



Andre et al. 



(120 lOf ). These authors report a typical filament 



and Polaris Flare regions by 
width of 271 ~ 10^ AU, or ~ 0.05 pc. In turn, filaments extending for over 3 pc can be 
readily seen in their Fig. 1, implying aspect ratios of up to ~ 60. From our Fig. [31 we 
see that this in turn implies filament c ollapse times up to ~ 10 times larger than the 



Andre et al. 



(I2OIOI ) report a typical column density 



corresponding spherical one. Moreover, 
of the filaments of ~ 10^^ cm~^. We can thus estimate the mean volume density in 
the filaments as n ~ N/2TZ ^ 6.5 x 10^ cm"^, for which, assuming a mean molecular 
mass fi = 2.36, the spherical iiee fall time is t^ sd ~ 1.3 x 10^ yr, and the Jeans length is 
Lj ^ 0.08 pc. According to our results, a filament of length 2L = 3 pc would thus take a 
time ~ 1.3 Myr to collapse. But if it contracts towards its center of mass, then quickly the 
central parts of the filament will reach a volume density large enough that a segment of the 
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filament of lengtfi equal to its tfiickness will become Jeans unstable. Indeed, the required 
volume density for the Jeans length to become Lj = 2TZ = 0.05 pc is nj ~ 1.8 x 10^ cm^^, 
or only ~ 2.8 times the mean volume density of the filament. At this point, the central 
clump becomes locally Jeans unstable, and proceeds to collapse on a timescale up to one 
order of magnitude shorter than the filament, naturally producing a star- forming core with 
a filamentary "appendix" that accretes onto it. Moreover, this will cause the core to sustain 
star formation for times significantly longer than its own (spherical) free-fall time, since its 
gas supply is continuously replenished by the accretion flow from the filament. 

In turn, filaments may form from the collapse of sheets. As mentioned in sec. [31 
during the collapse of a thin sheet, the acceleration is maximal at the edge, and thus 
the collapse occurs from the outside in, forming a dense ring in the periphery, which is 



topologically equivalent to a filament. Such hie rarchical fragmentation 



filanients has been suggested by several authors ( 



1984 



Hanawa et al. 



1993 



Kofman fc Poggsvan 



numerical simulations fjTurner et al. 



1995 



Schneider &: Elmegreen 



Tom sheets to 



1979 



Gaida et al 



1995), and has a. 



Burkert fc Hartmann 



so been readily observed in 



2004; 



Hartmann fc Burkert 



2007: 



Vazguez-Semadeni et al. 



Prieto et al. 



2007 



Heitsch et al. 



2009 



Rosas-Guevara et al 



201 



201l[ ). Thus, the collapse of giant molecular clouds may proceed in a way that 



is anything but a monolithic, three-dimensional collapse. 



5.2. Implications for the free-fall estimate of the SFR 



A second implication of our results is that they may contribute, at least partially, 
towards alleviating the well known "star formation conundrum". The latter consists in 



that the observed star formation rate (SFR), of a few Mq yr ^ (e.g. 



Smith et al. 



1978 



Diehl et al. 



2006 



Robitaille fc Whitnev 



2010l ). is roughly two orders of magnitude lower 



than the simple estimate obtained by dividing the total molecular mass in the Galactic ISM 
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by the three-dimensio nal free-fall time correspondi ng; to the mean density and temperature 



of the molecular gas (IZuckerman fc Palmer 



of the original proposal by 



19741) ■ T his argument caused the dismiss 



Goldreich &: Kwanl ( 119741 ) that molecular clouds should be 



in gravitational collapse, and was replaced by the notion that t 



correspond to supersonic microturbulence (IZuckerman fc Evans 



le ob served linewidths 
1974). The latter was 



thought to provide support against the clouds' self-gravity, a notion that persists until 
the present. Later on, support by magnetic fields was considered as well, and led to the 
notion that molecular clouds are in near virial equilibrium, with SF proceeding at a much 



slower rate than free- 
reviews bv IShu et al 



all, only as allowed by m ediation of ambipolar diffusion (see, e.g., the 



1987 



Vazquez-Semadeni et al. 



Mouschovias 



20001 : 



19911) and of lo cal turbulent compression (e.g.. 



Mac Low fc Klessen 



200J). 



However, if molecular clouds and their substructure are indeed in near free-fall, as sug 



gested by various recent studies 


[e.g., 


Burkert & Hartmann 


2004; 


VazQuez-Semadeni et al. 




2007. 


2010; 


Hartmann & Burkert 


2007; 


Peretto. Hennebelle & Andre 


2007; 


Galvan-Madrid et al. 


2009; 


Schneider et al. 


2010; 


Csen 


geri et al. 


201lJ 


b 


), then the SF conundrum reappears as a 



chal 



enge to theoretical models of free-falling clouds (e.g., 



Zamora-Aviles fc Vazquez-Semadeni 



20111 ). Our results may imply that, if the cold (molecular and atomic) gas is distributed in 



sheets and filaments rather than in three-dimensional structures, as suggested by various 



observational studie s (e.g., 



Molinari et al 



2010 



Ballv et al. 



Andre et al 



1989 



de Geus et al 



1990 



Heiles &: Troland 



2003 



20101), then the 3D value of the free-fall time is an 
underestimate to the real value, implying that the 3D free-fall SFR value overestimates its 
actual value, and that the severity of the SF conundrum may be reduced by up to one order 
of magnitude. 

Of course, the correction provided by our results is not expected to provide a full 
solution to the SFR conundrum, since it is unlikely that sheets and filaments have the 



- 16 - 



required aspect ratio values of A ~ a few xlO^ to completely account for a two-order- 
of-magnitude discrepancy between the observed and the free-fall estimate values of the 
SFR (see Figs. [T]and[2]). Besides, the presence of even weak magnetic fields (implying 
supercritical mass-to-magnetic flux molecular clouds) and stellar feedback are expected 
to reduce the SFR beyond whatever is achieved by the possible low dimensionality of the 
clouds. Nevertheless, our result may provide a reduction of the discrepancy factor between 
the observed SFR and its free-fall estimate by factors ranging from a few to almost one 
order of magnitude, since observed aspect ratios of sheets and filamen ts in the ISM can 



reach values of up to 10^ Meiles fc Troland 



2003 



Molinari et al 



20101 ) 



6. Summary and Conclusions 

We have computed the free fall time for two geometries that depart from spherical 
symmetry, namely, thin circular sheets and long filaments. In both cases, we considered 
two different behaviors for the volume density, one assuming that it remains constant 
through the collapse and the other assuming that it is the mass that remains constant, and 
the density increasing as the object collapses. These two cases should bracket the actual 
collapse timescale. 

We have parameterized the problem by the aspect ratio A of the object, defined either 
as the ratio of the cloud's radius to its (small) thickness {A = R/h) in the case of sheets, 
or as the ratio of the cloud's half-length to its (small) radius {A = L/TZ), in the case 
of filaments. Our calculations assume that the clouds contract only along their largest 
dimension (radially for sheets, and longitudinally for filaments), and thus our results are 
more accurate as larger values of A are considered. The value A = 1 is a frontier value 
between the two symmeties, although it is also the value for which our calculations are 
in largest error, since three-dimensional contraction should ensue in that case, with the 
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free-fall time being given by Tfr^sD- Nevertheless, even in this extreme case, our calculations 
only deviate from the 3D value by factors of order unity. 

For both cases, we found that the collapse time increases as A^l"^ . In particular, this 
imphes that the collapse time of a thin sheet is of the order of that of a sphere with much 
lower volume density, so that it is its column density that is the same as that of the sheet's. 
This result has two important implications for the structure of molecular clouds and star 
formation. First, it naturally explains the common morphology observed in molecular 
clouds, where star-forming or pre-stellar clumps are embedded within filaments that appear 
to be accreting onto them. According to our results, this is a natural consequence of the 
longer free-fall time for a filament than for any clump-like structure within it that contains 
enough mass to be itself collapsing. Because the filament has a longer collapse time, it 
will continue to accrete onto the locally collapsing 3D clump after the latter has managed 
to increase its density by a large enough amount as to become distinguishable from the 
filament. 

Second, it may imply that, if molecular gas in the Galaxy is distributed in primarily 
low-dimensional structures such as sheets and filaments, then the so-called SFR conundrum 
may not be as strong as it is normally stated, because the relevant free-fall time for the 
cold gas in the Galaxy may be longer than has been considered. However, we expect that 
this effect is hkely to only account for a fraction of the discrepancy between the observed 
and the free-fall prediction for the SFR, since the required aspect ratios to fully account for 
the conundrum would be too large (^4 ~ a few xlO^), and besides other effects are known 
to contribute to reduce the SFR, such as magnetic fields and stelUar feedback. In any case, 
our results suggest that determining the topology of molecular clouds is important for 
estimating their true expected collapse timescales. 

We acknowledge useful comments from an anonymous referee, which helped improving 
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CONACYT 102488 to E.V.-S., UNAM-DGAPA grant PAPIIT IN106511 to G.C.G, and a 
fellowship from CONACYT-SNI to J.A.T. 
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Fig. 3. — Dependence of the free-fall time for a filamentary cloud on the parameter A. The 
asymptotic behavior oc \/A is seen to be reached quickly as A is increased. 
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